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The field of very high energy (VHE) ground-based 7-ray astronomy, pioneered by the 



Whipple 10m atmospheric Cherenkov telescope (ACT) [I], is being revolutionized by the 
second generation observatories, H.E.S.S. [[2], MAGIC [ 3], and VERITAS [S]. Dozens of 
new 7-ray sources, some representing new classes of VHE emitters, have been discovered 
during the last decade [[5]. In space based gamma-ray astronomy a leap in technology from 
\q \ SAS-2/COS-B [El [7| to EGRET [[S] lead to an order of magnitude increase in detected 
sources; the leap to GLAST/LAT [EJ [TU] technology is anticipated to produce a further 
order of magnitude increase. The lesson to be learned is that increases in sensitivity, 
angular resolution and ultimately in scientific output comes through progress in detector 
technology. In anticipation of such progress it is valuable to re-evaluate the design limita- 
tions of current generation ACT observatories and to investigate where improvements can 
be made that will lead to the next generation. In this paper we focus on one component, 
the optical system. 

1. Introduction 

ACTs operate in a fundamentally photon starved regime. To avoid contamination of 
the Cherenkov image with night sky background photons, the exposure must be matched 
to the intrinsic duration of the Cherenkov light pulse, ~ 6 nsec for an atmospheric cas- 
cade. Therefore, unlike in conventional optical astronomy, the image cannot be improved 
through increased exposure. This motivates the development of optical systems with very 
large primary mirrors, having diameters in the range of 10 — 30 m. In addition, the optical 
systems of ACTs must be composed of the minimal number of optical elements, typically 
one, to circumvent light loss [ [111 121 131 E] . Thus, to image a Cherenkov light flash with 
a few nsec exposure, the optical system has to be "ultra-fast" , extending the terminology 
of photography. 

The light detector of an ACT records the individual photons produced as Cherenkov 
radiation by high energy secondary particles, which are themselves created as a result of 
the interaction of VHE 7-ray photons and cosmic rays (CRs) with the atmosphere. The 
atmospheric cascade, which is the physical source of light, develops at the altitude of < 15 
km, and the cascade core has a mean impact distance of a few hundred meters or more 
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with respect to the telescope position. Therefore, the imaging of Cherenkov radiation from 
a cascade demands a telescope imaging sensor which accommodates field angles of at least 
1.75°. There are several technological motivations to increase the field of view of future 
ACTs substantially above the minimum of 3.5° — 4°. A larger field of view improves the 
collecting area of ACTs for the highest energy 7-rays (E 7 > 10 TeV), since their cascades 
can be detected at larger field angles when they impact at a distance of several kilometers [ 
[To] . Observations of extended sources, such as supernova shells, galactic molecular clouds, 
or the hot plasma of galaxy clusters, in which 7-rays are produced through interactions 
of diffuse cosmic rays with ambient matter, would also benefit from increased field of 
view. At energies below a few hundred GeV the sensitivity of ACTs is limited by the high 
rate of atmospheric cascades caused by CRs impacting the atmosphere. To accurately 
estimate this background, and minimize systematic error, it is desirable to contain within 
the same field of view both the putative source and a few equivalent regions of empty 
sky. A field of view of ~ 10° or larger would accommodate observations of galactic 
and extragalactic extended sources. Finally, the detection of transient VHE phenomena, 
such as GRBs or rapid AGN flares, provides strong motivation to pursue a future ACT 
observatory capable of monitoring the full sky with high sensitivity. Currently, observing 
such transients requires an external trigger from a wide field of view survey instrument, 
almost exclusively an X-ray satellite with relatively small collecting area. In principle, a 
large array of ACTs, each having a field of view of 10° — 15°, would be able to provide 
coverage of ~ 1 — 2 steradians of the sky and generate a self-trigger to initiate more 
sensitive pointed observations [ITS]. 

Amongst optical telescopes, ACTs are characterized by an unsurpassed value of etendue 
(throughput), the volume of the phase space of incoming rays focused by the telescope, 
equal to the telescope primary mirror area times the field of view solid angle; its typical 
value for the present day ACTs is 2 x 10 3 m 2 deg 2 . The UK Schmidt Telescope (UKST) has 
the largest etendue among existing optical telescopes, ~ 60 m 2 deg 2 , only 50% larger than 
that of the Palomar Schmidt, ~ 40 m 2 deg 2 , built in 1948. The Large Synoptic Survey 
Telescope (LSST), a state of the art project under consideration, will have etendue of 
~ 250 m 2 deg 2 [ [TT] . A future ACT with 12 m diameter primary mirror and 15° field of 
view would exceed the etendue of LSST by a factor of 80. In view of this comparison, 
the only reason that ACTs remain relatively inexpensive instruments (the cost of LSST 
exceeds $100M) is that these telescopes have low angular resolution; several minutes of 
arc versus the 0.5 seconds of arc of LSST. 

Angular resolution has not been the main parameter driving the design of current 
generation ACTs and consequently the resolution of their image sensors has been limited 
to ~ 10 minutes of arc per pixel. At energies larger than a few TeV the sensitivity of ACTs 
depends weakly on angular resolution, giving little motivation for improvement of this 
parameter. However, at lower energies (< 1 TeV), the telescope operates in a background 
dominated regime, and the angular resolution has a major impact on its sensitivity. It 
has become evident from recent simulation studies that the angular resolution of current 
generation ACTs is far from the ultimate, dictated by the physics of Cherenkov cascades 
in the atmosphere [HH [16]. To accurately reconstruct the arrival direction of a 7-ray, the 
core of the atmospheric cascade, created by the secondary electrons and positrons with 
energies roughly above 100 MeV, must be identified in the shower image. Simulations 
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show that the transverse angular size of the core is ~ 1 minute of arc. Figure [T] illustrates 
that the reconstruction of the arrival direction of 7-rays deteriorates slowly when the pixel 
size of the imaging detector increases from 1 to 2 arc minutes, and more rapidly when 
the pixel size exceeds 4 minutes of arc. Thus, the optics of future ACTs should facilitate 
imaging with a pixel size of 1 — 4 arc minutes, a factor of 2 — 9 better than present day 
instruments such as VERITAS and H.E.S.S. (9 arcmin). For telescopes consisting of a 
single mirror surface, this requirement translates to optical systems with /-ratio (ratio of 
focal length of to primary mirror diameter) in the range 2.5 — 3.0, so that the aberrations 
are contained within the size of a pixel in the imaging detector. Such optical systems are 
significantly "slower" than those of current ACTs for which the /-ratio is 1 — 1.2. 

An optical system with both large aperture and large /-ratio leads unavoidably to a 
large camera plate scale. For example, a 10 m aperture telescope with / = 2.5 has plate 
scale of 0.73 cm per arcmin, and the desired resolution of a few arcmin matches well to the 
size of a few cm diameter photo-multiplier tube (PMT). Existing ACT technology utilizes 
expensive arrays of ~ 10 3 individual PMTs read out by complex, fast electronics with 
nsec timing. To improve the image resolution by a factor of ~ 3 and increase the field of 
view by a similar factor, the number of pixels must increase by two orders of magnitude. 
Issues of cost, reliability of operation, and the difficulty of maintaining 10 5 PMTs may 
limit the scalability of existing technology to future wide field of view instrumentation. It 
appears that highly integrated, reliable, and inexpensive imaging sensors will be required 
to address these problems. The discussion of possible technological solutions, such as 
multi-anode PMTs, Silicon PMs, image intensifiers, etc., is outside the scope of this 
paper. However, an important feature of all of these novel photon detecting devices is 
that their small physical size makes them incompatible with the large plate scale of single 
mirror optical designs. 

In this paper we explore the capabilities of two mirror optical systems for ACTs in 
an attempt to simultaneously balance the competing requirements of small plate scale, 
moderate tolerance to aberrations, large field of view and large telescope aperture. We 
set the scale of the optical system by the choice of the physical size of the pixel in the 
imaging detector, thereby fixing the physical dimensions of all elements, such as the 
diameter and focal lengths of the primary and secondary mirrors. This approach permits 
a direct comparison of the performance of single- and two-mirror optical systems assuming 
that the same imaging detector is installed in the focal plane of each. Rescaling of the 
pixel size proportionally changes the dimensions of all optical elements without affecting 
the conclusions of this study regarding the quality of imaging and the achievable field of 
view. In addition, we limit the amount of acceptable aberration to the angular size of 
the pixel, which has been constrained by the physics of atmospheric cascades. Although 
the smallest structures of the Cherenkov light images corresponds to a pixel size of ~ 1 
minute of arc, in this study we consider only 3' pixels giving an improvement of the image 
resolution by a factor of ~ 3 over present day ACTs. Without this compromise, the 
number of channels in the camera would exceed one million. With the parameters of the 
camera fixed we investigate one- and two-mirror optical systems to determine the largest 
effective light collecting area, given the desirable telescope field of view. 

The prescribed plate scale fully determines the focal length of the optical system, there- 
fore we are constrained to consider systems for which the product of the optical system 
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Figure 1. Ratio of the number of 7-ray events reconstructed within angular distance 9 of 
the source to the square root of the number of isotropic background events in the same 
region. The two figures show the results of simulations for 40 and 100 GeV photons. A 
hypothetical, future generation array of 75 m 2 telescopes with separation of 80 m and 
with image sensor quantum efficiency of 25% is assumed [H5J. The simulations employ 
the CORSIKA air-shower package [HI]. To isolate the effect of image sensor pixellation, 
an idealized telescope model, consisting of two reflective losses with zero dispersion due to 
PSF, was assumed. The curves show the effects that increasing pixellation in the camera 
have on the ability of the array to reconstruct the arrival direction of the primary 7-ray. 
For 1 arc minute pixels the signal-to-background ratio reaches a maximum at approxi- 
mately 7 and 3.8 minutes of arc from the source location for 40 and 100 GeV photons 
respectively. With 8 arc minute pixels, slightly smaller than those of current generation 
instruments such as VERITAS, the ratio reaches a maximum value at approximately 
10 and 5 minutes of arc. The relative amplitude of the peaks indicate reduction of the 
instrument sensitivity by ~ 40%. 
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/—ratio and the primary mirror diameter is constant, fD p = F = const. An increase 
of the aperture diameter to improve the light-gathering power unavoidably results in a 
corresponding decrease in the /—ratio, in turn amplifying all primary aberrations, such 
as spherical aberration ~ l// 3 , coma ~ 5/f 2 , as well as astigmatism and field curvature 
~ 5 2 /f. We define the "limiting aperture" as the maximal size of the primary mirror 
for which the optical system produces a high quality image of a point source located at 
the maximal field angle 5 max - For such an optical system, twice the RMS of the point 
spread function in the tangential or sagittal plane (whichever is the largest), 2a (<5 ma x), is 
equal to the angular size of a pixel in the camera. With a single-mirror optical system 
it is possible to correct spherical aberrations by using a continuous, parabolic primary 
or by utilizing the Davies-Cotton design, which has a discontinuous primary mirror. In 
these systems coma limits the maximal aperture of the telescope. With the two-mirror 
system, both spherical aberrations and coma can be corrected (Schwarzschild theorem [ 
enabling a further increase of the telescope aperture. One such aplanatic (free of 
spherical aberrations and coma) Cassegrain system, known as Ritchey-Chretien [ 20J, is a 
frequent choice for optical telescopes designed to achieve a high image quality over a large 
field of view, e.g. the Keck telescopes, VLT, the Hubble Space Telescope, WIYN, etc. 
In optical telescope applications the secondary mirror is configured to magnify the image 
and increase the optical system focal length. In contrast, ACTs call for a reduction of 
the image by the secondary in order to achieve a decreased focal length. The traditional 
Ritchey-Chretien configuration of optical telescopes utilizes hyperbolic mirrors for both 
the primary and secondary to correct first order coma over field angles of a few minutes of 
arc. The correction of coma over a larger field of view, as required for ACT applications, 
was investigated by Chretien in his original works on the Ritchey-Chretien design [ |2"T]. 
We follow his approach, considering optical systems which correct for coma over a field 
of view of several degrees by departing from hyperbolic primary and secondary mirror 
surfaces. In this paper we evaluate the benefits of such an aplanatic telescope (A-T) for 
ground-based 7-ray astronomy. 

The paper is structured as follows. In section [2] we discuss the most traditional single- 
mirror ACT design, the Davies-Cotton, to establish a benchmark for further comparisons. 
Section [3] is devoted to the explanation of the A-T optical system. In section H] we provide 
the results of ray-tracing simulations performed to determine the optimal parameters for 
a two-mirror aplanatic telescope. Section [5] contains detailed characteristics of the optical 
system optimized for application to ground-based gamma-ray astronomy and provides 
a summary of the results. Throughout this paper we deliberately avoid discussing the 
advantages and disadvantages of the D-C and A-T optical systems from the perspective 
of design, fabrication and alignment of the mechanical system and of cost. The latter 
depends strongly on the specific scientific goals which the instrument must achieve and 
on the changing costs of different components of the telescope. In the discussion, section [HI 
we briefly address both subjects. 

2. Ideal Davies-Cotton telescope 

To illustrate the deficiencies of single-mirror telescopes for wide-field imaging applica- 
tions we review the properties of the Davies-Cotton reflector, which is the most commonly 
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used design in ground-based 7-ray observatories. Originally, the Davies-Cotton telescope 
was developed as a solar concentrator [Ell, an d as such, it does not satisfy the rigorous 
requirements of astronomy in the visible wavelength range. Nevertheless, the design has 
been widely appreciated by 7-ray astronomers for three primary reasons. A large reflector 
composed of many small, identical, spherical facets is relatively inexpensive to build. The 
alignment of the optical system is trivial. The performance at large field angles is better 
than that of a single spherical or parabolic reflector [[TJ. A Davies-Cotton telescope con- 
sists of a spherical primary mirror, f(tp,6), with radius equal to the focal length of the 
system, 

(sin 9 cos tp 
sin 9 sin tp 
1 -cos# 

The normals of the individual facets, n, are aligned to a single point (x, y, z) = (0, 0, 2F) 
so that 



n 



sm 2 cos tp 
sin 2 sin tp 



COS | 




Since the normals to the facets do not coincide with the normals of the surface of the 
reflector, the ideal Davies-Cotton design is discontinuous at every point of the primary 
mirror and therefore the facet size of an ideal telescope must be equal to zero. The effects 
of finite facet size, which have been investigated through detailed ray tracing [[251 126] . 
degrade the imaging by introducing an additional source of aberrations: the astigmatism 
of individual facets. We omit these effects, and consider only global aberrations, which are 
irreducible even in ideal Davies-Cotton optics. 

For the purposes of ray-tracing we denote the direction of an incoming photon as 
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where 5 is the field angle. The direction of the ray reflected from the primary is given by, 

(— cos 5 sin 9 cos tp + (l — (1 — cos6>) (cos tp) 2 ) sin S 
— cos S sin 9 sin tp — (1 — cos 9) sin if cos tp sin S 
cos 5 cos 9 + sin 9 cos cp sin 5 

The crossing of the focal plane by the reflected ray, 

sin 9 cos p \ 
sin 9 sin tp + ts, 
4 Ff ^ / \ 1 - cos 9 1 

requires that the parameter t be given by, 

cosfl 5 2 F 1 / r /sin0\ \ ^, r4 , 

cos o cos 9 + sm 9 cos tp sm o 4 i*/ cos 9 \ \ cos 9 J J 
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To assess the effect of the field curvature we introduced the i"f^ 2 term, assuming that 
2Ff is the radius of the curvature of the focal plane. The formulas are valid to the third 
order in the field angle. The coordinates of the intersection with the focal plane are given 
by: 

* • , 1 + (r^b - 1) cosV 

— = sino- 



F cos 5 + ^| cos ip sin 5 

cos ~ 

+ T~^~~77 (- sin # cos y? + <5 (l + (— - - l) cos 2 cp 
4 Ff cos \ \ \ cos J 

V_ = ^ t^b" 1 ) (sinv?) (C0SV?) 
F cos 5 + cos <z> sin 5 

COS y " 

5 2 F 1 / / 1 \ 

+ "T ^ a ( ~ sin 9 sin V + 6 ( q - 1 ( sin <f) ( cos V 5 ) 

4 F/ cos \ \ cos J 

The relevant moments of the light distribution in the focal plane of the telescope are 
determined by averaging over the reflector aperture, 



Jo 



'(27) I* / o du) 
/ (...)8/ 2 sm#cos^-^-, 

JO 27T 



where the primary /-ratio, / = |J, is assumed to be much larger than 1/2 (expansion 
parameter is jp). The position of the image centroid in tangential (x-axis) and sagittal 
(y-axis) coordinates is given by: 

7) = °- 



5 4 ( ( F\ 2 if 49 F F 2 
+ ^ 4+— +— 35 + —— + -, 



The centered second moments of the light distribution are equal to: 

2 l 0/ 4 ^ 4/ 2 

2^(( 4+ ^>; ' 6/ 2 r" ' AF f ' F| 



4) 2 



+ 




>|)(a1)) = o 



8 



V. Vassiliev , S. Fegan and P. Brousseau 



The contributions from the five primary aberrations can be read from these formulas 
characterizing the distribution of light on the focal plane. The spherical aberration term, 
oc l// 6 , is eliminated in the Davies-Cotton reflector design. The coma term, oc S 2 /f A , 
is the dominant source of aberrations for all field angles 5 < -j. For a maximal field 

angle, 5 max expressed in degrees, this implies that / < f j*— (/ < 2 for <5 max > 7.16° 

and / < 3 for <5 max > 4.77°). Coma cannot be corrected by introducing curvature into 
the image plane. Detailed investigations of various single mirror prime-focus designs for 
wide field of view IACTs [ 27] suggest that the Davies-Cotton reflector is among those 
with the smallest coma. For example, for a continuous parabolic primary the tangential 
RMS due to coma is a factor of \[2 larger, \/2^j S , while the sagittal RMS is similar to a 
Davies-Cotton reflector. A continuous spherical primary would add spherical aberrations 
and increase tangential coma by a factor a/14/3 = 2.16. A hyperbolic primary reduces 
tangential coma, however the sagittal component remains unaffected. Astigmatism, the 
next aberration term, oc 5 A j f 2 (RMS oc 6 2 /f), has the same scaling with 5 and / as 
the effect of field curvature. Either tangential or sagittal astigmatism can be corrected 
by the proper choice of curved focal plane. For instance, the choice of infinite radius 
of curvature, 2Ff — > oo, eliminates sagittal astigmatism, and the choice of 2Ff = — \F 
cancels tangential astigmatism (the minus indicates that the focal plane is curved towards 
the primary). If the focal plane has curvature 2Ff = —F, the overall blurring caused by 
astigmatism is minimized but not completely eliminated. The curvature of the focal plane 
can also affect distortion, the fifth source of aberration, which arises from the change of 
the plate scale with field angle (<5 3 term in (jlA formula). Although the distortion can be 
dramatically reduced by the choice 2Ff = — ^F, in practice the reduction of astigmatism 
is preferred since distortion can be accounted for through calibration. 

We consider wide field of view Davies-Cotton telescopes in which coma is the dominant 
source of aberration. To contain the image of a point source at 5 max within a single pixel 
of the camera sensor, of size p minutes of arc, one needs 




2 4 / 2 



5„ 



1° 



P 



1 

< 60 Ll'J 



which implies / > | y [npp-] j ■ The diameter of the camera, d = [^pp] jgfD, exceeds 
half of the primary mirror diameter (i.e. > 25% obstruction of mirror area) when £ max > 

2 1 

5° [|^] 3 [|t] 3 . Thus, the Davis-Cotton design appears to be an option for a telescope 
with wide field of view, not exceeding 10° — 11°, for pixel size p of ~ 3 minutes of arc. A 
slightly smaller field of view (by ~ 10%) can also be achieved with the parabolic mirror [ 
[2Tj . However, these designs require an optical system with large /-ratio (> 2.5), for which 
the physical size of the pixels in the detector, = -^f [£] must exceed ~ 2 cm for a 
10 m mirror aperture. This large plate scale prohibits the use of small, highly integrated 
image sensors, and appears to leave only the option of using an array of ~ 4 x 10 4 

2 r 'i 2 

[ 5°^] p" individual PMTs. More cost effective and, perhaps, robust technological 
solutions, such as multi-anode PMTs with a typical plate scale of 3 — 6 mm per pixel, or 
even more integrated light sensors such as SiPMs, CMOS, etc., cannot be accommodated 
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in this design without additional optical elements, such as lenses and light concentrators. 
This limitation and the desire to break the "10°" field of view limit (<5 max = 5°) motivates 
our study of two-mirror telescope designs which may address both issues simultaneously 
Catadioptric systems, with a large aperture lens in front of the camera, are disfavored for 
ACT application because the technique depends on the detection of broadband Cherenkov 
radiation, most of which is concentrated in the blue and UV regions of the spectrum. 
These systems would suffer from considerable chromatic aberrations, and be subject to 
significant absorption around 300 nm, which is a characteristic of acrylic or glass Fresnel 
lenses. 

With fixed plate scale, — , the maximal primary mirror radius for which the comatic 
aberrations at <5 max field angle are contained within the size of a single pixel is given by 



D 
~2 



lm 



3.6 



7T 



P 



5mm 



5° 



(2) 



The effective mirror area, including obstruction of the primary by the camera, is 



.4 



(Sn 



DCeff 10 

max 



4.125m 2 
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-1 2 



5mm 



5° 
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"3'~ 
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18 




V 





5° 



(3) 



This analytical result neglects the change of the PSF (equations [Q) due to obscuration. 
However, since the PSF is dominated by rays reflected from the outer regions of the mirror, 
obscuration of the inner portions by the camera is not expected to modify our conclusion 
considerably. Equation [3] provides a benchmark to which different telescope configurations 
with angular pixel size p, coupled to physical pixels of size P, can be compared. 

3. Aplanatic two-mirror telescopes 

3.1. Definition of optical elements 

The basic geometry of the optical system is illustrated in figure [2j Throughout this 
section we introduce a number of variables to describe aplanatic two-mirror telescopes. 
The most important are summarized in table [TJ 



3.1.1. Primary Mirror 

The surface of the primary mirror is given by 



r(ip, h) 



hcosip 
h sin (p 
z(h) 



(4) 



where <p G [0, 2tt] is the polar angle, h G 



^ 



is the running surface radius, D„ is the 



primary mirror diameter, and z (h) is the sag function, that can be expanded as a Taylor 
series 



F p Z 1 



h 



+ Z 2 



h 



h 



2d 



+ ... 



(5) 
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Figure 2. Schematic diagram showing the A-T design and the parameters used to describe 
the focal lengths and diameters, and placement of the optical elements. 
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Table 1 

Review of variables used in the text. 



Variable 


Eqn. 


Description 


7 ~» 7— 1 

V> P 




Diameter and local length oi primary mirror, bee hgure |2J 


£ 

fp 




i-ratio oi the primary f p = r p /D p . 




4 


Polar angle on the surface oi the primary. 


7 

h 


4 


Running radius on the surface of the primary, h G [U, ±J p /2\. 


z(h) 


4 


Sag function of the primary. 


Zh 


6 


Derivative of primary sag function, Zh = dz/dh. 


D s , F s 




Diameter and focal length of secondary mirror. See figure [zj 


£ 

Is 




f-ratio of the secondary f s = F s /D s . 


J. 


n 1 1 

ecu 


rolar angle on the surface of the secondary. 


P 


rm 


Running radius on the surface of the secondary, p G [0, D s /2]. 




rm 

LUJ 






fTHl 


Derivative of secondarv sat? function ti = du/dn 


a 


EH 


Position of the secondary with respect to the primary, in units 






of the primary focal length. See figure [2j 






Diameter and focal plane curvature. See figure [2j 


v> 


da 


Polar angle on the focal plane. 


A 


us 


Running radius on the surface of the focal plane, A G [0,Df/2]. 


/(A) 


da 


Field curvature function describing the focal plane. 


P,p 




Physical and angular size of a pixel in the camera. 


F 




Focal length of the combined optical system. 


I+77 


7 


De-magnification due to secondary, relating F and F p . 






Field angle, angle of incoming rays with respect to optic axis. 


<^max 




Maximum field angle, defining field-of-view. 


5 ~~ a(l-H?) 


E3 


Critical parameter of the solution for primary sag function. 


X 2/p 


[1011221 


Critical parameter of A-T, related to the maximal opening angle 






of the rays impacting the focal plane. 
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By convention we choose Z\ — \ so that F p is the focal length of the primary mirror. 
The dimensionless argument ■#■ runs from zero to =4-, where f p = jf-is primary mirror 



2/p' 

/—ratio. If the primary were chosen to be a conic surface, then the constants Z n would 
be determined as expansion coefficients of 



(1+« P ) „2 
4 J/ 



1 + K» 



In this case, the conic constant k p is related to the coefficient Z 2 as k p = 2 6 Z 2 — 1, and the 
remaining aspherical coefficients are determined completely. The normal to the primary 
surface is 



(6) 




dz ( h 

Zh= dh =2Zl [K 



h 



h 



2n-l 



Assuming the final focal length of the two-mirror telescope is F, we parametrize the 
magnification of the secondary mirror by 



1 



l + ?7 



(7) 



so that the desired image reduction is achieved when rj > 0. Since the plate scale is fixed 
by the given physical and angular pixel size, P and p respectively, the primary mirror 
diameter must satisfy 



V 



7T 



60 180 



D r 



l + r? 



7T 



3.6 x 10 3 



P 



5mm 



5 x 10~ 3 m, 



and therefore 



5m 



3.6 



7T 



5mm 



1+7/ 

2/ P 



(8) 



In this definition 77 and f p are free, unit-less parameters of the design and the physical 
dimensions of all optical elements of the A-T system can be re-scaled by the appropriate 
choice of P. In analogy with the Davies-Cotton reflector (equation [3]) the effective area 
of the optical system can be expressed as 



.4 



ATeflf 10 max 



103.1m 2 



p 


2 


"3 /_ 


2 




5mm 




.P. 


( 


V 2/ p J 



Tr (f p , 77, S r 



(9) 



where Tr is the light transmission coefficient, accounting for obstruction of the primary 
mirror by the secondary and the secondary mirror by the camera, and for vignetting due 
to the loss of some oblique rays entering optical system at 5 max which miss the secondary 
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mirror. Section @] describes the results of numerical simulations of ^ATeff (^max) for various 
design parameters. 
The ratio, 

appears frequently in the following calculations and, for this reason, we identify its approx- 
imate physical meaning. According to Liouville's theorem etendue is preserved through 
the optical system. Requiring that the etendue at the primary mirror and at the focal 
plane be equal, 

(D p V /<5 max 71 V ( F p 5 max 7T \ 2 / . $ X2 

7T — - X 7T ~ 7T X 7T 2 Sm 



2 J V 10 180 / \1 + V 10 180 / V 2 

suggests that this ratio is related to the maximal opening angle, 2$, of the rays impacting 
focal plane x = Yf 1 = 2 sin ^. 

3.1.2. Secondary mirror 

The surface of the secondary mirror is defined analogously 

/ p cos (f) \ / \ 
f(0,p) = psin0 + , (11) 



PCOS0 ^ 




S 


psin0 


l + l 




v(p) J 







where p is the running surface radius, is the polar angle, v (p) is the sag, and aF p is 
the primary-secondary separation. Following the traditional Cassegrain design of two- 
mirror systems, the secondary mirror is placed inside of the primary mirror focus (a < 1). 
However unlike the Cassegrain design, which uses a magnifying, convex secondary mirror, 
we use a concave reflector to reduce the size of the image. The optically equivalent, 
Gregorian-like configuration (a > 1) with convex secondary mirror is disfavored because 
it leads to a larger primary-secondary separation and ultimately to a more expensive 
telescope. Reduction of the image by the secondary mirror, by the factor 

F 1 F 

_ = 1 = 1£ (12) 

F p I + 77 F,-{l-a)F p ' 1 ; 

is achieved with the use of concave secondary with negative focal length, 

F, = Jl^Ht, (13) 
V 

indicating that it is curved toward the primary. The sag function v (p) can be expanded 

as 

-'■(«fe) ,+ «fe) i+ -*fer-)- <i4) 

where V± — \. Again, if the secondary is chosen to be a conic surface, its conic constant 
k s would be given by n s = I^V^ — 1. The normal to the secondary mirror is expressed as 

v p cos (f) 

T , v p s\n<P I (15) 

V 1 + v p 
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where 



v 



dv 
dp 



2Vi 



+ W 2 



+ ...2nV n 



2n-l 



+ 



An approximation of the diameter of the secondary mirror is given by 



D ps 



D r . 



r-t \ , n O max 

l-oc)+af p — — 



The exact value of D s can deviate significantly from this estimate. The diameter of the 
secondary mirror for zero field (5 max — 0) is uniquely determined by the configuration of 
the optical system, through the parameters a and D p . The actual diameter of a practical 
system, however, may exceed this value and should be determined through optimization 
in which the effects of obscuration, vignetting, and aberrations are balanced. We discuss 
this subject further in section [51 



3.1.3. Focal Plane 

The distance to the focal plane from the secondary mirror 
optical ray tracing. With respect to the primary mirror, the surface of the focal plane is 
determined by 



F p follows from on-axis 



r(ip, A) 



/ A cos ip ^ 




( 







I A sin ip 


I* 








\ /(*) J 




\ 







(16) 



where A is the running surface radius, ip is the polar angle, and / (A) is the field curvature 
function, which can be expressed as a Taylor expansion in the form 



/ 




A 



+ Y 2 — +...Y n — 



A 



2d 



+ 



(17) 



Again, the parameter Y\ = 1/4 is chosen by convention. The radius of field curvature, 
2Ff, and, if necessary, the high order aspherical terms Y 2 = ^ (1 + «/), etc. should be 
chosen to minimize astigmatism. In section H] we describe a numerical optimization which 
finds the tangential and sagittal focal planes as well as the optimal focal plane, which 
minimizes the maximal contribution to astigmatism from both projections. 

3.2. Constraints on the optical system 

Assuming that the photon beam is incident on the primary mirror at an angle S in the 
xz— plane, 

sin S 

'1= | o 

— cos 5 



the direction of the ray reflected at r(ip, h) is given by 

-2z h cos^\ sin5 ( l-z£cos(2y>) 
-2^siny3 | + _ 2 , - | -z^sin(2y9) 

2zh cos <f 



9 ~ 2n(<f, ft) 



cos 5 
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The intersection of a ray reflected from the primary with the secondary, 



/ hcosip ^ 






h sin cp 




[( 


V * ) 







— 2z ft COS<p 

—2z h simp I cos 5 + 



1 - z 2 cos (2p) 
-z\ sin (2y>) 
2^ cos <p 







^ p COS \ 


| sin 5 


H 


p sin 






\v + otF p J 



requires that parameter t be given by 

aF p + -u — z 
(1 — 2^) cos 5 + 2zh (cos (p) sin 5 

For an on-axis ray (5 = 0) the solution necessitates that (p 



and 



p = h — (otFp + v 



>2 ' 



(19) 



The ray, reflected from the secondary mirror, propagates to the focus of the telescope 
along a path described by the unit vector 



1 



H + 1) (zi 



1 



2 (f p - Z h ) (1 + 2^) cos ip 
2 (u p - z h ) (1 + ZfcUp) simp 

(f/s - - (1 + Z h V p f 



Therefore an optical system which is free from spherical aberrations must satisfy the 
following requirement, 



1 (l-q) F 



2 (v p - z h ) (1 + z h v p ) 

{Vp - Z h f - (1 + Zftfp)" 



(20) 



The mathematical conditions for correction of coma can be derived from ray tracing 
through a perturbative Taylor expansion of equation [IS] with respect to 5. The decep- 
tively simple "Abbe sine condition", that ensures the exact cancellation of the comatic 
aberration term in the expansion, requires that the optical system have a particular kind 
of symmetry. Namely, the intersection points of all rays reaching the focus of the optical 
system with the incoming rays which produced them (beam parallel to the optical axis) 
must form spherical surface, i.e. 



P 



h 



U ^ 1 + 7} r P 



VR 2 - h 2 

According to equation [19] in the limit of small h 

v 2z h 



(21) 



p = h — (aFp + v 



h 



1 



h - (aF p ) AZi-r + O (h 2 ) = (l-a)h + (h 2 ) 



the radius of the spherical surface can be identified with the focal length of the two-mirror 

system, 
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To satisfy the Abbe sine condition for all running h from to D p /2 one needs D p /2 < 
F P I (1 + 77). For a coma free optical system, the reduction of the image by the secondary 
mirror is therefore limited to 

This inequality prohibits the construction of two-mirror systems with arbitrarily small 
plate scale or alternatively arbitrarily large aperture. As a direct consequence, the tele- 
scope imaging camera diameter cannot be made smaller than D p ^^^-. 

The surfaces of both the primary and secondary mirrors of an ideal aplanatic optical 
system are completely constrained by equations [191, [201 an d ED If we denote 

G = . 1 andVG 2 - 1 = (1 + 77) ^-G 

'l -(<! + ->>*)' 



then 



_ {l-zl)h-2z h ((a-^f)F p - 



p = J(P . 

2z h + (1 - z\) VG^l 

= (1_- z\) h (1 - (1 - o) GO - 2z h (aF p - z) 



2z h + (1 - 4) Vg^T 



By substituting for the variable h and re-scaling the sag functions and p one obtains, 

h 7 (1,7, 

Vy=a + v) Yp ,Z=(l + V )-, z h = - = 2^Z y , (23) 

(i + v) _ ^ (I - yg) - ^ ((" (i + g) - (i - gO) - g) (0A s 

~P — yV ry n I (\ r/i\ ' \ ' 



;i + 77) _ (J-^J)(vT^l/-(l-a))-Z I ,(a(l + 7 ? )-Z)vT^i/ 



1.' 



^ 3Vr=y + (J - y^j) 

(i(TTv^) + ( 1 + - /r ^)^-^) 



(25) 



p 



VV 7 n n 7W\ • ( 26 ) 

ZyVi -y- \\-y z l) 



These equations can be reduced to a single second order non-linear differential equation 
for the primary mirror sag function, 

z _ 1 {{a{2 + V )-Z-l)Z y + yZl-\) 

yy 2 (a(2 + r ] )-Z-l + VT^) (l-y + JT=y)' 
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which must be solved with the initial conditions Z\, n = and Z v \ n = \ T , \ T . We found 

ly=0 y\y=0 4 (I+77) 

that the solution, Z (y), has a unique, non-trivial property, which is likely related to the 
presence of a hidden symmetry in aplanatic optical systems, allowing it to be written 



a 



Z=-qy-(l-a)V(y,q), 



(27) 



where g _1 = a (1 + 77) and ^ (y, q) is a function of only two independent parameters 
rather than the expected three: y, a, and (1 +rj). This degeneracy requires that be 
a solution of a first order differential equation which can be readily obtained after some 
mathematical transformations 



dy 



The solution satisfying ^\ y=0 = and expressed in quadratures is given by 



* (y, q) 



1 

1 

(2- 



1 (g(y/r^-l)+2) 1+ ^ 

2 (v^i/ + ir 1+I ^(>/r^i/-i) : 



/ 



q) (s 



(s + l) 2 ~^(g (s-l) + 2) 2 ^ i-« 



or 



2 (VT^i+l)- 1+ - V 2 / (,(VT^-l) + 2) 2+ ^ 
A non-divergent solution for \l/ (y, g) exists if 

2 - = 2a(l + r ] )>l-VT^= X " 

5 1 + V1-X 2 



■dx .(28) 



or 



1 + v 713 ^ 



(29) 



Analytical solutions can be found for q — 1 — — , where m is any integer. Simulations 
indicate that the optimal designs are always found for q in the range of [|, |] and for this 
reason two boundary solutions have particular importance 



*(y, g = l/2) 



32 



*(j/,g = 2/3) = 



162 



(v^i+l) : 



(l9(l-y) + 28v^i+7) 
(v^^+1) 2 



(30) 



(31) 
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For the range of the numerical optimization discussed in the following section, the solution 
(y, q — 2/3) always provides an effective telescope mirror area which is within a few 

percent from the maximal achievable value. 

Within its region of convergence, the general series solution can be represented in the 

form, 



^3 
\1>4 

^5 
\[>7 



( 4 + <?) 
gi i g(30 + llg + 2g 2 ) 

j^mq(336 + U6q + 4lq 2 + Qq 3 ) 

r q (5040 + 2414g + 829g 2 + 194g 3 + 24g 4 ) 



(32) 



2949,120 
82 575 360 



q (95 040 + 48 504g + 18 754g 2 + 5489g 3 + 1114g 4 + 120g 5 ) . 



Since 

Zi = 
Z n = 



CM 



:n > 1 



it is evident that the primary mirror is not a conic surface for any choice of a and 77. 
For the Cassegrain-like design (a < 1) the primary is, to first order, a concave hyperbolic 
mirror in with conic constant K p given by (1 + K p ) = —2 ( ^ — l) (1 + rf) 2 . However, for all 
values of a and r/ it deviates from a hyperbolic surface in the second expansion order. For 
the Gregorian-like configuration (a > 1) the deviation of the primary from an elliptical 
surface is also unavoidable for any physical parameters a > 0, and rj + 1 > 0. 

The parametric definition of the secondary mirror can be computed numerically based 
on the equations [231 EH and [25j We have not been able to derive a general analytic 
solution representing the mirror surface. The series solution, assuming it is convergent, 
is given by: 



M£) +V2 U) +v 3 U) +v 4 U) + ... 



ikv-Hv + i) 3 (12^ + 1(^-11) + ^ 

/ 11 o . 



^3 = ^ 



^=«k*r 7 fa + i) 4 ( 



v $ = vdmv 9 fa + K 



1 20(^-1(147/ + ^ -77) 



/ 



(577 - 3) + ^ 
1680 7 



V 



+| (369r/ + 49r/ 2 + 3r/ 3 - 1357) 
(778r/ + 29r/ 2 + 189) 
+i (6I77 - 19) + £ ' 



\ 



/ 



(33) 



V S = sgltilo^ 11 fa + 



/ - 30240 (^TT) 

(17487/ + 327t/ 2 + 38t/ 3 + 2r/ 4 - 4843) 

"+§ (3258r/ + 219r/ 2 + 5r/ 3 + 1604) 

'* +4(131r] 2 - 590r/-181) 

V %^(29^-7) + ^ 



/ 
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The secondary mirror is also not a conic surface. It can be made to approximate a 
conic surface in the first and second order by a specific choice of a and 77 such that 

{h 2 + 777 + 5) (2r] + 3) 

and a — 



rf + § ?7 2 " " (77 + 1) 77' 



For a G [§, l] this requires that \ (v / 13 + l) < V < \ + 3) and -1.704 < k s < 
—0.742. Numerical optimization suggest that the preferred values of 77 are in the range 
[1,2.5], outside of this bound. Thus, like the primary, an optimized secondary mirror 
deviates from a conic surface in the second order. 



4. Numerical ray-tracing simulations 

The basic imaging characteristics for various optical system (OS) configurations were 
obtained through ray-tracing simulations. Each simulated OS has the same plate scale, 
5 mm per 3 arc minutes, and therefore the same focal length of 5.73 m, and is characterized 
by three parameters a, f p , and x- The shape of the primary mirror was calculated by 
numerical integration of equation [28] and the surface of the secondary mirror was derived 
through the parametric representation, equations [231 EH arid [251 This method finds the 
diameter of the secondary mirror under the zero field approximation. The parameters were 
varied within the following ranges a G [0.45, 0.85], f p G [0.7, 3.2], and \ £ [0-60, 0.98]. For 
each configuration the curvature and the conic constant of the focal plane were numerically 
computed by minimizing astigmatism. For this purpose, the surfaces of minimal tangential 
and sagittal astigmatism were determined and the focal plane was defined so that the 
effective diameter of the PSF, A ps f = 2 x max{RMS sag i t tai, RMS ta ngentiai} is minimized. 
In the majority of configurations tested the tangential RMS dominates over most of the 
range of field angles and hence the focal plane is effectively determined by minimizing 
tangential astigmatism. It is typical that the tangential RMS exceeds the sagittal by 
10-15% for field angles larger than 3 degress. For smaller angles, the saggital RMS may 
exceed the tangential slightly, but they are always within a few percent of each other. 

With the focal plane resolved, the effective PSF diameter, A ps f, and the effective light 
collecting area, A e s, were computed as a function of field angle, 5, in the range 5 G [0°, 7°]. 
In these calculations, both obscuration by the secondary mirror and by the 14° FoV camera 
were accounted for. The configuration with the largest A cS was determined, subject to 
the condition that A ps f at various field angles does not exceed 5, 4, 3, and 2 minutes 
of arc. It was found that all optimal configurations have af p ~ 1 (af p G [0.8,1.2]). 
For these systems the primary to secondary separation is within 20% of the diameter of 
the primary mirror. Optical systems with smaller separations are disfavored due to their 
large obscuration ratio, D s /D p , and those with larger separations suffer from considerable 
vignetting, since the secondary mirror is not large enough to intercept all oblique rays 
reflected by the primary. We observed that the performance of the optical systems was not 
a rapidly varying function of the defining parameters a, f p , and \. Therefore the bounds 
on the volume in this three dimensional space which contains the optimal configurations 
are not strict. Our goal is to narrow the parameter space which needs to be scanned to 
find optimal solutions. If for the purposes of engineering and implementation, an optical 
system must deviate slightly from the best configurations, the bounds can be relaxed. 
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A non-trivial finding is that the parameter q = [a (1 + f])}~ , which was of prime impor- 
tance in the derivation of the optical surfaces, is constrained to a narrow interval q G [5, §1 
for all favored solutions. For the majority of these systems q is closer to the upper bound, 
therefore the surface of the primary mirror can be described approximately by the ana- 
lytic solution, equations [27] and [3TJ Finally, for all optimal configurations the parameter 
X was found to be within the range x [0.70,0.97]. The larger x values correspond to 
optical systems with larger primary mirrors, larger A^, but significantly increased A ps f. 
In the opposite limit, the point spread function is considerably improved, at the price of 
decreased effective light collecting area. While analyzing the optimal solutions, a corre- 
lation between the x an d I + rj parameters was observed, suggesting that their product 
remains nearly constant, 2.2 < %(1 +77) < 2.5. This relationship implies that systems 
with a larger degree of de-magnification at the secondary mirror must have larger f p , and 
consequently smaller primary mirror diameter and smaller opening angle of rays at the 
focal plane. 

The performance of the A-T and D-C designs is compared in figure |3j The left panel 
shows A,ff(<5 max ) which can be achieved by the designs, if they have the same plate scale 
and A ps f(5 max ) is limited to 3 minutes of arc for both. The D-C curve is given by equa- 
tion El the A-T curve is compiled from the characteristics of the best designs (a, f p , x) 
obtained through simulation. It is evident that the decreased aberrations of the aplanatic 
telescope allows for an increase of A eS (S max ) by a factor larger than 10. 

In practice, two effects will reduce the effective light gathering power of the A-T design 
compared with the single mirror D-C design. First, a fraction of the light will be lost at the 
second reflection. Current generation ACTs have mirror reflectivity in excess of 80% at 350 
nm, even under conditions of extreme weathering. The industrially coated and protected 
mirrors used in astronomical facilities routinely have reflectivity in excess of 95%. Second, 
the larger angles at which the photons impact the light detecting sensor in the A-T design, 
up to 45° — 50°, may introduce light loss. Utilization of anti-reflection technology, such as 
a thin film of semiconductor deposited on the entrance window, can significantly reduce 
the reflection coefficient. The large opening angle of the rays impacting the focal plane in 
an aplanatic telescope may also reduce the efficiency of light concentrators placed in front 
of the photon sensor. However, the fraction of light lost likely does not exceed 10 — 20%. 
We discuss this subject further in section |6j These factors combined do not come close 
to the order of magnitude (or larger) increase in the effective light gathering power of the 
A-T system. 

The effect of vignetting was studied as a final stage of numerical optimization. To 
extend the diameter of the secondary mirror beyond the zero field approximation, the 
series solution given by equation [321 and [331 was used. It was found that for the majority of 
investigated configurations, a series solution for the primary and secondary sag functions 
with six terms provides a sufficiently accurate description of the mirror surfaces such that 
Ap S f (5) remains unaffected. However, the Taylor series for the optical systems with the 
largest A e g (6) converge relatively slowly and the sixth-order series approximation was 
determined insufficient to reproduce the PSF. For this reason, the numerical solutions 
obtained directly from the differential equations were fit with sixth-order polynomials and 
these fits were then used to study vignetting. Increasing the diameter of the secondary 
mirror reduces A eS (0). For oblique rays, 5^0, the total change of A ef[ (5) is due to 
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the combination of two competing effects. Although increasing the obscuration ratio 
reduces A e g (S) , the additional rays intercepted by the secondary mirror compensate for 
this reduction at sufficiently large 5. Figure @] shows that the effective light collecting area 
can be equalized over a large range of field angles. However, increasing the diameter of 
the secondary degrades the PSF of the telescope and therefore limits the range of field 
angles for which a uniform response can be achieved. 

Ultimately, the particular configuration appropriate for any application must be de- 
termined by balancing the requirements for uniformity of A e R (5) over a given range of 
5, maximizing effective area, A e ff (0), and minimizing PSF, A ps f (5). Table [2] shows three 
configurations tailored to different goals. OS 1 is optimized to achieve the largest effective 
area, with relaxed requirements for the compensation of vignetting and imaging quality. 
OS 2 is designed to provide the largest field of view with a uniform response, which is 
accomplished through a significant reduction in effective light gathering power. For this 
system the secondary mirror is extended by 25% and the inner region of the primary 
mirror area is removed. For the design of OS 3, we required both a field of view 14 de- 
grees and an effective PSF diameter of less than three minutes of arc. The performance 
characteristics of all three systems are shown in figure [5j 

5. Telescope performance characteristics 

A cross-sectional view of OS 2 showing the primary, secondary and focal surfaces is 
given in figure This figure also illustrates ray tracing for field angles of zero and five 
degrees. OS 2 provides a compromise between a relatively large effective light gathering 
area and a large field of view. This configuration is also unique for the degree to which 
it allows for compensation of vignetting. The coefficients of the primary and secondary 
mirrors for this configuration are given by 

Zi = {0.25, -0.189377, -0.604706, -4.21374, 21.8275, -425.160} 
V t = {0.25, 0.013625, -0.010453, 0.014241, -0.012213, 0.005184}. 

Both the primary and secondary mirrors can be segmented to reduce the cost of the optical 
system. A possible arrangement of mirror facets, as "petals" , is shown in figure UJ This 
scheme has the advantage of requiring a minimal number of different surface shapes. A 
study of the tolerance of alignment and positioning of mirrors is beyond the scope of this 
paper. Nevertheless, our experience with the simulations suggests that the requirements 
are stricter than those applied to the H.E.S.S. and VERITAS optical systems. The use 
of automated alignment and calibration systems will likely be required, e.g. [|2"9"]. 

A composite showing the light distribution at different field angles of OS 2 is given in 
figure [H The structure in the distribution is shown in more detail for two field angles in 
figure [9j A considerable distortion is evident in the first figure. As a result of minimizing 
the astigmatism, the plate scale becomes a slowly varying function of field angle, which 
is expressed as a cubic term in the mean of the light distribution on the focal plane, as 
illustrated in section [2J In astronomical telescopes the effect of distortion is removed 
during the processing of images. For ACTs a similar approach remains valid, although 
additional care will be required to account for a slightly variable density of night sky 
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Figure 3. Left: The maximal effective light collecting area of a D-C and an aplanatic 
telescope, each constrained to have a plate scale of 5 mm per 3 minutes of arc and 
A ps f(<5 max ) < 3'. Right: The ratio of A e s(5 max ) for aplanatic and D-C telescopes. The 
three curves correspond to A ps f(<5 max ) constrained to be less than 3, 4, and 5 minutes of 
arc. 




Field angle [deg] Field angle [deg] 



Figure 4. Left: The effective area as a function of field angle for OS 2, table [2j The 
dashed line shows A eS (5) when the secondary has been extended by 25% to compensate 
for vignetting. Right: The effective PSF diameter, A psf (5), as a function of field angle. 
The extension of the secondary mirror of OS 2 successfully equalizes the effective area 
for field angles up to five degrees with almost no effect on the PSF. This property is 
not generic for all OSs; normally the secondary mirror can be increased by only ~ 10% 
without significant degradation of imaging. 
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Table 2 

The parameters of three optical systems, tailored for the largest effective area (OS 1), a 
uniform response over the largest range of field angles (OS 2), and the minimal aberration 
(OS 3). 



Parameters 


AC 1 


AC 9 




n 

Up 


1 n on m 
lU.zu m 


O /ID m 

y.4U m 


O 09 ™ 

o.oz m 


n 

JSp inner 


A D7 m 
^.Ul 111 


^.UO 111 


A 91 m 

'i.ZO 111 


Primary Mirror Area 


68.70 m 2 


52.15 m 2 


47.04 m 2 


fp 


1.46 


1.80 


2.08 


F P 


14.892 m 


16.915 m 


18. 346 m 


a 


0.66 


0.59 


0.56 


l + r] 


2.599 


2.952 


3.2032 


Secondary Position 


9.829 m 


9.980 m 


10. 274 m 


D s (Zerofield) 


5.35 m 


5.29 m 


5.03 m 


D s 


5.77 m 


6.61 m 


6.23 m 


D s inner 


0.60 m 


0.50 m 


0.40 m 


Secondary Mirror Area 


25.87 m 2 


34.12 m 2 


30.36 m 2 


F s 


-3.167 m 


—3. 553 m 


-3.664 m 


Focal Plane Position 


7.881 m 


7. 631 m 


7.754 m 


Ff 


-1.313 m 


-1.481 m 


-1.479 m 


K f 


-1.00 


-4.28 


-6.62 


$ = 2 arcsin 


52.8° 


48.4° 


45.3° 


Plate Scale 


10 cm per deg 


10 cm per deg 


10 cm per deg 


F 


5.73 m 


5.73 m 


5.73 m 


A e fT (Equalized) 


55.26 m 2 


34.92 m 2 


30.54 m 2 


A cS (Vignetted, Zerofield) 


59.85 m 2 


47.61 m 2 


41.24 m 2 


Apsf (5°) 


5.60' 


3.32' 


2.33' 
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Field angle [deg] Field angle [deg] 



Figure 5. Left: The effective function of field angle for the three configurations 

of OS summarized in table [21 Right: The effective diameter of the PSF of the light 
distribution in the focal plane of OSs. 



5n 1 1 1 1 1 1 1 1 1 1 1 n 5^ 




Cross-sectional plane Z [m] Cross-sectional plane Z [m] 



Figure 6. Illustration of incoming rays traced through the optical system to the focal 
plane for tangential rays at field angles of zero (left) and five (right) degrees. 




Figure 7. One possible scheme for faceting the primary and secondary mirrors. Four 
different facet types are used on the primary, three on the secondary. Each facet is 
limited to an area less than approximately 0.45 m 2 with no linear dimension larger than 
1 m. 




2 3 4 5 
Focal plane X [deg] 



8 



Figure 8. Composite of images made in the focal plane of an aplanatic telescope with 
configuration OS 2 by rays at field angles between zero and seven degrees. 
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background photons in the field of view of the telescope. An analogous effect will also be 
introduced if vignetting is not compensated. 

Another characteristic important for ACT applications is the time dispersion introduced 
by the telescope, which must be less than the intrinsic spread of arrival times of Cherenkov 
photons (a few nanoseconds). For a D-C telescope, such as VERITAS or H.E.S.S., a planar 
front incident on the telescope results in an almost uniform spread of the photons at the 
focal plane, over the interval of 4 — 5 nanoseconds. Such telescopes are not isochronous, 
even for on-axis photons. The time dispersion increases further for oblique rays and 
scales linearly with the diameter of the primary mirror. In comparison, the A-T design is 
isochronous for on-axis rays, and only a small time dispersion is introduced when photons 
enter the OS at large field angles. Figure [10] shows the results of a simulation of this 
effect. 

6. Discussion 

We have presented a study of the design of a new, wide-field, two-mirror, aplanatic tele- 
scope for ground-based 7-ray astronomy. It has been shown that this design provides a 
potential alternative to prime-focus optical systems utilized by present-day observatories. 
If the same plate scale is used in both optical systems, the A-T design radically outper- 
forms the D-C in terms of effective light gathering power, ability to accommodate wide 
field of view, and amount of time dispersion introduced by the telescope to the arrival 
times of Cherenkov photons. 

Up to field angles of five degrees, both designs offer solution for an optical system with 
equal light collecting area and an effective PSF diameter smaller than 3 minutes of arc. 
Although this paper does not present a detailed analysis of the advantages and disadvan- 
tages of the mechanical properties of the D-C and A-T designs, some are immediately 
evident. The D-C design requires a very large focal ratio and photon detector plate scale. 
Assuming an effective light collecting area of 25 m 2 (D p = 6.2 m), field of view of 10° 
{Smax = 5°), and moderate angular pixel size of 5 minutes of arc, one finds that the re- 
quired focal length is 11.5 m (f p = 1.88). This translates to a camera of diameter of 2 
m with 10 4 pixels, each of which must be a 1.8 cm diameter PMT. Such a camera would 
inevitably have a large mass, and require a sturdy support structure, which would itself 
present a challenging engineering problem. The lowest fundamental frequencies of the 11.5 
m length camera support structure would need to be suppressed to avoid degradation of 
the optical performance. The moment of inertia of the telescope will likely be dominated 
by the heavy camera, increasing the time to slew and stabilize the system. However, 
the clear advantage of the D-C design is the simplicity of the fabrication, assembly and 
alignment of the reflector. 

The A-T design can be made more compact. The equivalent effective light collecting 
area of 25 m 2 can be achieved with a primary mirror of 7.5 m diameter, after compensating 
for obscuration effects and loss of light due to the second reflection. The secondary mirror 
of 4.25 m diameter must be supported at the distance of 7.2 m from the primary. The A-T 
design allows for rigid support of the camera through the 3 m diameter hole in the primary 
mirror. The significantly smaller plate scale of the A-T design allows packing of the 10 4 
6 mm pixels within a 75 cm diameter camera. The lowest natural frequencies of the A-T 
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Figure 9. Images made in the focal plane of aplanatic telescope with OS 2 configuration 
by rays at field angles of zero (left) and five (right) degrees. The origin is chosen at the 
image centroid. 
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Figure 10. Left: RMS spread in arrival time of rays at the focal plane of aplanatic 
telescope with OS 2 configuration as a function of field angle. Right: Distribution of 
arrival times of photons impacting the OS at a field angle of 5 degrees. 
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optical system are almost twice that of the equivalent D-C configuration. The mass of the 
primary and secondary mirrors is distributed almost on a spherical surface minimizing 
moment of inertia. This may allow for fast slewing drive systems and alternative cross- 
elevation axis mounts. The A-T design allows for a field of view wider than ~ 10 degrees, 
with pixel size less than 3 minutes of arc, while any prime-focus optical system will suffer 
from an excessively large obscuration of the primary mirror by the camera. However, 
the clear disadvantage of A-T is the complexity of the optical system consisting of non- 
spherical, off-axis mirror elements, challenging mechanical support of the large secondary 
mirror, and much tighter tolerance requirements for primary and secondary support and 
alignment. The use of a mirror support structure in which the deflections of the primary 
and secondary mirrors are matched to maintain correct alignment is likely necessary to 
overcome the problem of telescope flexure. Such a structure might use the Serrurier truss 
design, used commonly since the 1930s to mount the primary and secondary mirrors in 
large optical telescopes, such as the 6.5m MMT and 10m Keck telescopes. The unavoidable 
need to reduce the OS weight and its moment of inertia will probably require the use of 
modern, light-weight, high-strength materials, such as carbon fiber-reinforced plastics, 
to construct the optical support structure. Alignment of the segmented optical system 
may require automated edge sensors and actuators, since the tolerances of the aplanatic 
telescope OS are roughly equivalent to a mm range radio telescope, such as ALMA [EB], in 
the range 100-20 /xm. Although these requirements are three to four orders of magnitude 
looser than the diffraction limit to which the positions of the segmented mirrors of the 
Keck telescopes are maintained, they are far stricter than those of current ACT optical 
systems, such as VERITAS, MAGIC, and H.E.S.S. Many of the required technological 
solutions are, however, already being explored by various groups, e.g. [[2H1 USB ED I3"2"] . 

In the traditional D-C design the focal plane is not curved, since the main source 
of aberration, coma, cannot be corrected by field curvature. The flat imaging camera is 
assembled from individual PMTs, allowing approximately 75% of the camera to be covered 
by photon sensitive area (e.g. the 19 mm diameter Hamamatsu R3479, with 15 mm 
diameter photocathode, arranged in a hexagonal pattern would satisfy the requirements 
of the D-C telescope described in the example). Light concentrators, typically Winston 
cones, are used in front of each PMT to collect additional photons from the dead space. 
The efficiency of light concentrators grows as the opening angle at which photons impact 
the focal plane is decreased. Collecting efficiencies of ~ 50% for dead space photons are 
routinely achieved by the present day ACTs; higher efficiencies may be possible in optical 
systems with larger f p . Thus for D-C telescopes, light loss at the focal plane can be as 
small as ~ 10%. It is apparent that light concentrators for use on aplanatic telescopes 
cannot be as efficient, since the full opening angle of rays at the focal plane is large, in the 
range 90 — 100 degrees. However, the plate scale of the A-T design is significantly smaller, 
6 mm per pixel in the example given, and this allows the use of a more highly integrated 
light detectors. As an example, the Hamamatsu H8500 8 x 8 MAPMT has a dead space 
of ~ 16%. Light concentrators, though less efficient, can still be used to recover some 
of the lost photons bringing the difference between the A-T and D-C optical systems 
to at most a few percent. In fact, small cones can be used in the A-T design to help 
define the curved focal plane and therefore minimize astigmatism. The use of SiPMs and 
HPDs, with potentially higher quantum efficiency but significantly larger dead space, may 
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present a problem. The curved focal plane of the A-T design can be coupled with large 
aperture electrostatic image intensifiers, since the focal plane is naturally convex, e.g. [ 
E21 Ell- I n this implementation the dead space is automatically zero, if the photocathode 
has a uniform quantum efficiency over the full entrance window of the intensifier. For this 
type of the photodetector the angular resolution of the aplanatic telescope at small field 
angles is compatible with an image sensor which has an angular pixel size smaller than 
one minute of arc in the central part of the camera, if desired. It is possible to install 
a flat camera in the focal plane of an A-T system at the cost of slight degradation of 
imaging due to increased astigmatism. 

Throughout this paper we have intentionally avoided discussing the relative costs of the 
D-C and A-T designs. Putting aside the fact that the cost of the individual components 
of ACTs are constantly changing, and that any estimate would be quickly out of date, 
a meaningful cost comparison can only be made when the instrument design has been 
decided upon through optimization to address particular scientific goals. For example, 
one could optimize the instrument to conduct a survey of the largest possible area of sky, 
at the photon energies around ITeV, in the shortest amount of time possible to achieve 
a predetermined sensitivity. Alternatively, one could optimize for the highest sensitivity 
to point or extended sources at energies above ~ 10 TeV, or to improve the performance 
in the lowest energy domain. Thus the question of comparison of cost has neither an 
immediately apparent nor a unique answer. 

As an illustration of general considerations of an instrument cost estimation, we examine 
one possible strategy which could be used to optimize the parameters of a large array of 
telescopes, for which the footprint of the array is much larger than the characteristic size 
of the Cherenkov light pool. We also require that the array operates in the background 
dominated regime. One possible optimization parameter is the total throughput of the 
observatory, which is the volume of the phase space equal to the effective collecting area 
for 7-rays times the solid angle of the field of view. With our assumptions, the effective 
collecting area coincides roughly with the observatory footprint, i.e. the product of the 
number of telescopes, N t , and some characteristic area per telescope, A , determined by 
the choice of the telescope spacing. The latter is usually selected to balance the desire to 
increase the baseline for stereoscopic observations, giving a better angular reconstruction 
of an event, and to decrease it in order to achieve a for higher 7-ray trigger rate, especially 
at the lowest energies of interest. Thus, the throughput is oc A x N t X 7r x <5^ ax . We 
also account for effects of improved angular reconstruction due to the decrease of the 
detector pixel size, p, in the optimization parameter. A natural modification is T = 
A x N t x 7r x 5^ ax x A 2 (p), where A(p) expresses the values of the peaks of curves shown 
in figure [T], relative to the case of p — 1 minute of arc. Optimization of T is equivalent 
to the optimization of the sensitivity of the array over its field of view, if the effect of 
increased pixel size on the efficiency of background rejection is neglected. An expression 
for the approximate cost of the observatory is $totai = N t x ($ scope + N p x ($ pix + $d aq )), 



number of detector pixels per telescope, is the cost of the detector per pixel and $d a q 
is the average cost of the data acquisition electronics per pixel. 

Let us first assume that the cost of the camera dominates the cost of the telescope, 



where $ 



scope 



is the cost of positioner, optical support 
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Pixel size, p [arc minute] 



Figure 11. Various expressions involving A(p), the function describing the degradation 
in the sensitivity of a large array with increased camera pixellation, p, calculated from 
figured] (100 GeV). The function A(p) is normalized with A(larcmin) = 1. The various 
curves are approximately proportional to the "throughput" of a hypothetical instrument, 
T(p), under different sets of assumptions, as described in the text. The curves obtained 
for 40 GeV photons show similar trends. 
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§sco P e <C N p x ($pi x + $daq)- Then, with N t x N p constant to keep the overal cost fixed, 
the modified throughput, T = A X N t x ~ x N p X p 2 A 2 (p), is a function of angular 
pixel size only. Figure [TT] shows the behavior of p 2 A 2 (p) and indicates that the most cost 
effective strategy is to degrade the angular pixel size to p > 16 arc minutes, reducing the 
number of pixels per camera and degrading the angular resolution of the instrument, while 
improving the sensitivity due to the increased number of telescopes in the array. This 
conclusion is a direct consequence of our assumptions, that the array be large compared 
with the Cherenkov light pool and that it must operate in the background dominated 
regime, which were in turn required by a specific science goal. A different conclusion 
follows if observatory footprint is much smaller than its 7-ray collecting area and the 
array consists of a relatively small number of telescopes. In this case the collecting area 

1 fl 

and T scale with ~ N t rather than N t , and the constraint on overall cost requires that 
pA 2 (p) be maximized. The curves shown in figure [TT] cannot be used for quantitative 
conclusions in this limit, since A 2 (p) was not simulated for an array consisting of a small 
number of telescopes. Nevertheless, the qualitative trend is evident; the optimal angular 
pixel size is likely constrained within the interval p G [4, 10] minutes of arc. 

In the opposite limit, $ scop e 3> N p x ($ pix + $d aq ), unless the cost of the optical system 
and positioner is a faster growing function of the PSF diameter, A ps f (£ max ) = P, than 
A 2 (p), the strategy to improve sensitivity through improved angular resolution is prefer- 
able. It appears then that under the given assumptions, a cost efficient telescope design 
requires, roughly, that the cost of the camera doesn't exceed the cost of the optical sys- 
tem and positioner. A design that satisfies this criterion and, at the same time, achieves 
the ultimate angular resolution of ~ 1 minute of arc might be considered as "perfect". 
Current ACTs have pixels approximately a factor 8-10 larger than the "perfect" design 
and technological solutions enabling significant reduction of cost per channel, $ P j x + $d aq , 
are of great interest. 

The A-T design, with significantly smaller plate scale compared to the D-C design of 
equivalent effective light gathering area, may offer such an opportunity. To illustrate this 
we use the example discussed above in which D-C and A-T designs require pixels of 1.8 
cm and 0.6 cm respectively. The D-C telescope requires that the camera be implemented 
as an array of PMTs (e.g. Hamamatsu R3479), while the aplanatic telescope can also 
utilize MAPMT technology (e.g. Hamamatsu H8500). At current prices the per pixel 
cost of a single PMT exceeds that of a MAPMT by a factor of ~ 13 [[35]. The cost of the 
data acquisition electronics for the A-T design is equal to or less than for the D-C design. 
Both designs will benefit from the successful application of fast multiplexing, either at the 
telescope trigger level, which enables only the fraction of the full camera which contains 
the air-shower image to be digitized and stored, or at the digitization level where optical 
fiber delay lines are utilized to multiplex signals from several pixels into single FADC 
channel [ 36J. These technologies potentially have the promise to reduce $d aq by a factor 
of ten or larger. Furthermore, due to the small size of the camera in the example given, 
75 cm (29.5 in, 10°) - 40.6 cm (16 in 5.4°), the A-T design might be compatible with a 
camera utilizing a large aperture electrostatic image intensifier and a specialized 512 x 512 
or 1024 x 1024 pixel CMOS sensor [[37]. Although this technology has not yet matured 
far enough to offer a practical solution, it promises a substantial reduction in camera cost 
and a drastic improvement in reliability of camera operation, especially if applied to an 
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observatory with a large number (> 100) of telescopes. 

Camera cost can also be decreased through a reduction in field of view, equivalent to 
a reduction in the number of pixels, N p , while keeping the angular pixel size constant. 
The throughput, T, is invariant if the field of view is decreased by a factor of two and 
the number of telescopes is increased by a factor of four. Assuming that the cost of the 
camera does not dominate the total cost of the telescope, the strategy of increasing the 
field of view is cost efficient if the optical system and positioner for a single wide field 
of view instrument can be made less expensive than four instruments of half the field 
of view. For a modest increase in field of view over the 5 degrees currently utilized by 
ACTs, it is likely that the D-C design will continue to provide a cost effective strategy. 
However, a substantial increase in the field of view, and hence in the f-ratio, plate scale 
and the physical size and weight of the camera, will have a large impact on the cost of a 
D-C instrument, and may eventually become the dominant factor in the estimation of the 
total cost of the array. For an aplanatic telescope, a substantial up-front investment will 
be required for the development of the novel optical system, with its different, aspherical, 
off-axis mirror elements, which must satisfy the surface shape and alignment tolerances 
much tighter than presently required for ACT applications. At the same time, the wide- 
field A-T design potentially offers considerable savings through a dramatically reduced 
camera cost, and the reduction it allows in the number of optical systems and positioners 
required to cover the required solid angle in the sky, compared to the narrower field of 
view D-C telescopes. Evidently, an inexpensive process to manufacture and assemble the 
segmented optics of the aplanatic telescope is critical to making this design cost effective. 
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